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We show that the ground state energy of the translationally invari- 
ant Nelson model, describing a particle coupled to a relativistic field of 
massless bosons, is an analytic function of the coupling constant and 
the total momentum. We derive an explicit expression for the ground 
state energy which is used to determine the effective mass. 

1 Introduction and Results 

We consider the so-called translationally invariant massless Nelson model 
in d dimensions, which describes a quantum mechanical particle interacting 
with a quantized field of relativistic massless bosons [27] . Both, particle and 
field move in d dimensions and we assume that d > 3. 

The Hamiltonian of this model will be denoted by H\, with A 6 R denot- 
ing the coupling constant. The generators of translations in R , which are 
also known as operator of total momentum, commute with the Hamiltonian. 
This yields a direct integral decomposition of the Hamiltonian 
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with fiber Hamiltonian H\(P), [32]. This paper is devoted to the properties 
of the ground state energy 



E X (P) :=mia(H x (P)) 



(1.2) 



as a function of the coupling constant A and the total momentum P. Strictly 
speaking, one should refer to (jl.2p as the minimal energy, since in three 
dimensions the fiber Hamiltonian does not have a ground state [15] (for 
non-relativistic QED see also [23] ) . On the basis of heuristic physical prin- 
ciples, E\(-) should be twice continuously differentiable in a neighborhood 
of the origin, with positive second derivative. Establishing these properties 
is technically difficult due to the fact that the ground state energy is at the 
bottom of the continuous spectrum. By rotation invariance E\(P) only de- 
pends on the absolute value of P. Using functional integrals one can show 
that E\(-) has a global minimum at the origin, [21] , Frohlich showed in his 
thesis that the gradient of E\(-) is Lipschitz continuous [ISCE], see also 
|28j . Using renormalization [TJ [12] or iterated perturbation theory [17] it 
was shown for non-relativistic QED, that E\(-) is in a neighborhood of the 
origin a C 2 -function with positive second derivative; the same methods and 
results should apply for the Nelson model as well. However, these results on 
the regularity do not go beyond C 2 . 

In this paper it is shown that E\(P) is in a neighborhood of zero a real 
analytic function of A as well as P. We provide explicit analytic expansions 
of the ground state energy in powers of A. From this expansion it follows 
that on any closed subinterval of (—1,1) and for sufficiently small values 
of the coupling constant, the ground state energy E\(-) is a C 2 -function 
with positive second derivative. Furthermore, we derive a convergent small- 
coupling expansion for the effective mass 



in powers of A (see [37] for the study of the effective mass for the polaron 
and [221 E31 E] for non-relativistic QED). As is well known, Nelson's model 
studied in this article is very close to the polaron model which is of con- 
siderable physical interest (see [13] for a recent review). The literature on 
the polaron is enormous and a considerable part of it is concerned with the 
calculation of this effective mass (see \18\ d] for reviews on this particular 
problem). Some of the references most relevant to our approach for the 
calculation of m c g (A) are [25j [TOj, [361 US [22]. In view of the importance of 
such calculations, we therefore present an explicit combinatorial convergent 
expansion for the effective mass, in the Nelson model case. 



m e ff(A) : 



1 



u = l, ...,d, 



(1.3) 
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The method used in this article is a variant of the one introduced in [2] 
for the treatment of the massless spin-boson model. Our method exploits 
the well known formula which relates the minimal energy to the vacuum 
expectation of the semigroup generated by the Hamiltonian, i.e., 

E X (P) = Km -i Iog(n, e~ TH ^n). (1.4) 

T— >oo 1 

Using functional integration and a forest interpolation formula due to 
Brydges, Kennedy, Abdesselam, and Rivasseau (BKAR formula) we analyze 
the right hand side of (jl.4p . The functional integral representation of the 
quantity (f2, e~ THx ( p '£l) essentially provides a mapping of the original model 
to a statistical mechanics model of interacting intervals on the real line. 
The right-hand side of (]1.4|) then becomes the infinite volume limit of the 
pressure for this gas of intervals. The latter is then studied using a cluster 
expansion obtained from the BKAR formula which decouples the two-body 
interactions of these intervals in a minimal way. Namely, the formula can be 
viewed as a resummation in terms of trees of a sum over graphs with edges 
corresponding to the two-body interactions between intervals. This yields 
an explicit expansion of the minimal energy. 

We believe that the analyticity part of our result could alternatively 
be obtained using techniques based on operator theoretic renormalization, 
[U El fl2| 120) . However, the explicit expansion which we derive seems to be 
specific to the statistical mechanics methods employed in this paper. 

To illustrate the methods used in this paper we shall first consider a 
harmonic oscillator coupled to a relativistic field of massless bosons. This 
can be viewed as a didactic introduction to the full-fledged expansion in 
the translation invariant case. It will be shown that also for that harmonic 
oscillator model the ground state energy is an analytic function of the cou- 
pling constant A and an explicit formula for the ground state energy will 
be derived. We note that the analyticity part of this result can in fact be 
deduced for a larger class of confining potentials, not just a harmonic oscil- 
lator potential, using the result of Griesemer and Hasler in [20]. That result 
is based on operator theoretic renormalization and requires the Hamiltonian 
to satisfy a mild IR-condition. By working in an IR-regular representation 
of the CCR, [BJ [30], the Hamiltonian has sufficiently regular IR-behavior, 
for [20] to be applicable, while its minimal energy remains unchanged (we 
thank the referee for pointing this out to us). 

It would be interesting to further investigate whether the methods and 
estimates used in this paper can be extended to the renormalized Nelson 
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model, i.e., where the ultraviolet cutoff is removed according to the proce- 
dure outlined in |27j . 

Below we introduce the model and state the results. The bosonic Fock 
space is given by 

oo 

with vacuum vector $7 = (1,0,0, ...) over the Hilbert space f) = L 2 (R d ). By 
a*(k) and a(k), k E M. d , we denote the usual bosonic creation and annihila- 
tion operators in T satisfying canonical commutation relations 

[a(k),a(k')] = [a*(k),a*(k')]=0, [a(k),a*(k')] = 6(k - k'), (1.5) 

and 

a(k)Q = 0. (1.6) 

Equations (jl.5p and (|1.6p are understood in the sense of distributions over 
M d . We shall consider massless bosons with dispersion relation uj{k) = \k\. 
The operator of the free field is given by the self-adjoint operator 



H f := / oj(k)a*{k)a(k)d a k. 

JR d 

We define the real Hilbert space H-\i% to be the completion of 5 rea i 
with respect to the inner product 



' J Rd 2u)(k) 

For / G H_i/2 we define the field operator 

r fjdi, ^ ^ 
<f>(f):= - 7 ==(f(k)a(k)+f(-k)a*(k)). 

It is well known, see for example [8], that the following estimate follows as 

1/2 

a consequence of the Cauchy-Schwarz inequality. For all <p G D(Hj ) 

< 2||//o;||(,||F} / V|| + ll/M^IMI, (1.7) 

where || • |L denotes the norm of f). 

Let p G i?_i/2 and p/w G f). For define p x (-) := p(- + x). The 

Hamiltonian of the translationally invariant Nelson model is given by 

H x ■= --A x + H f + \(i>(p x ), 
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and acts in the Hilbert space H := L 2 (M. d )(£>J-. The Hamiltonian is invariant 
under translations and thus commutes with —iV x + Pt the operator of total 
momentum, where 

P f ■= [ ka*(k)a(k)d d k. 

JR d 

This yields the direct integral decomposition (jl.ip with fiber Hamiltonian 
H x (P) = ^(P-P f ) 2 + \ ( j ) (p )+H f , 

acting in T . By fjl -Tj) one sees that 4>(p x ) is infinitesimally small with respect 
to Hf. Hence H\ and H\(P) are self-adjoint on the natural domain of Hq 
and Hq(P), respectively [31] , We now state the main result of this paper. 

Theorem 1.1. Assume that p > a.e.. Then E\(P) is a real analytic 
function of A and P on the set 

|A|<ix< 1H P|)-lx(/ M A^f)". (1.8) 

The expansion coefficients are given in Equation (14.12H . 

Proof. First assume that p > a.e.. In that case the conclusion of Theorem 
11.11 follows as a consequence of Theorems 12.21 and 13-H Proposition 14.21 and 
Inequality (|4.17p . To extend it to the case p > 0, we choose an approximat- 
ing sequence (p n )n£N in with p n > f) n +\ > 0, such that for s = 1, 2 
the following limit converges in f), 

^ _£ n— >oo ^ _£ 
p n 0J 2 )■ puj 2 . 

Let E njX (P) := inf a(H (P) + \<t>{pn))- Then by ^HD) it follows that 

lim E n , x (P) = E X (P). 

n— >oo 

On the other hand, in view of ()4.17j) we can interchange the limit n — > oo 
with the infinite summation in the expansion (I4.12p of E nt \(P). Now each 
coefficient in that expansion converges as n — > oo, by dominated conver- 
gence. □ 

Remark 1.2. The assumption p > in Theorem \l.l\ includes the case which 
is of physical interest. Nevertheless, one could relax this assumption for 
example as follows. Let Xo-(^) = l[o-,oo)(I^D- The assertion of Theorem \1.1\ 
holds, provided for all A and P satisfying (j 1 .8 j) the Hamiltonian ^(P—Pf) 2 + 
\(f>([Xa * pj) + Hf has for any a > a ground state which has a nonzero 
overlap with the vacuum. 
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Let us now consider a particle in a harmonic oscillator potential coupled 
to the quantized field. The Hamiltonian of this model acts in rl and is given 
by 

L x := H osc + Hf + \<t>{p x ), 

where 

_1 1 2 _rf 

-"osc ; — 2 2' E 2' 

By (II. 7h . La is self-adjoint on the natural domain of Lo- We note that the 
operator is different from the explicitly solvable operator, which has been 
investigated in [5]. The next theorem states that the minimal energy of L\ 
is an analytic function of A. The region for which we can prove analyticity 
will depend on the quantity 

A := sup ( [ d d k \p{k)\ 2 \k\ n ~ 2 \ ' . (1.9) 

n>l \7]R d / 

Theorem 1.3. The infimum of the spectrum E\ := inf o~{L\) is a real 
analytic function of A for 



2s — i 
' Z > 2 



|A| < (2eA 

The expansion coefficients are given in Equation (|4,6 



Theorem 11.31 will follow as a consequence of Theorem 12.31 Theorem 
and Subsection 14.11 



2 Positivity Preserving Representations 

In this section we justify (|1.4|> . which will be the content of Theorem 12.21 
and we justify an analogous formula for the harmonic oscillator potential, 
which will be the content of Theorem 12.31 The following lemma will be used 
in the proof of Theorem 12.21 and Lemma 12.51 

Lemma 2.1. Let H be a self-adjoint operator on a Hilbert space 7~L. For 
ip € % we have 

lim -— log(V>, e~ TH ip) = inf supp ^ 

T— >oo 1 

where jjL^ denotes the spectral measure of tp with respect to the operator H. 
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Proof. Let = inf supp ^ . By the spectral theorem 

POD 

(V>,e" T ^)= / e- TX d^{\) 
and hence for any e > 0, 

where C e = fi^,([E^, E^ + e]). Taking the logarithm, we obtain 

^ < -- logty, e" TH V) < + e - - log C £ . 
The Lemma now follows after taking the limit T —> co. □ 
Theorem 2.2. Assume that p > a.e. and |P| < 1, then (jl.4p holds. 

Proof. W.l.o.g. assume that A > (for A = notice that f2 is the ground 
state of Hq(P)). Let T > 0. We write fj = f) rca i © if) rea ij where f) rca i = 
L 2 (M a! ;]R). We introduce ^(flreai) the real Fock space over f) rca i. This is a 
real Hilbert space. We define a Hilbert cone [H] in J-"(f) r eai) by 

oo 

C := 0CW := {/ G bf e iil(-l) re / > 0}. 

n=0 

It was shown in [151 [To] that e~ THx ^ p > is ergodic with respect to C (see also 
Section 3.3 in [33] and Propositions 2 and 3 in [13] )• It follows from the proof 
of Theorem XIII. 44 (c)=^(e) in [32] that e~ THx ^ is positivity improving. 
Fix e > and choose a nonzero function / with 

< / < e- eHx{p) fl, (2.1) 

where the inequality is understood with respect to the cone C, [13] • Now 
using that e~ THx ^ is positivity preserving we find 

-liogCn.e-P'+^Wn) < -Ii og( /, e -^)/). 

Taking the limit as T — > oo we obtain using Lemma 12.11 

inf a(H\(P)) < inf supp (J,q < inf supp jjLf. 

Let X denote the set of functions / satisfying (|2.ip . Since e~ eHx ^ P ' is pos- 
itivity improving, it follows that the linear span of X is dense in J 7 (I)). 
Hence inf o~(H\(P)) = inf inf supp fj,f, and the theorem follows using 
Lemma l2.ll □ 
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In the remaining part of this section we consider the harmonic oscillator 
potential and prove the following theorem. Recall that the ground state of 

H osc is (fo(x) := 7r~4e~^". 
Theorem 2.3. The following holds 

E x := mfa(Lx) = lim -- log(p ® e" TiA ^o ® fi). (2.2) 

To show Theorem 12.31 we introduce the Schrodinger representation of 
Fock space. For the following construction of measure spaces we refer the 
reader to [34] , By Minlos' theorem, there exists a measure d[i on Q := 
5^ eal (M d ) and a Gaussian random process £(/) on Q indexed by / G H_i/ 2 
with mean zero and covariance 

M(/X(s)) = (/.ff)-l/2. 

Moreover, there exists a unitary transformation f/jr : J 7 — >■ L 2 (Q,dfi) with 

t^n = i, u^(f)u^ = £(/). 

Define 

C/ p : L 2 QR d ) -> L 2 (M d , ^z), / ^ f/ Vo , 
which is a unitary transformation of Hilbert spaces. The map 

U:=U P ®U T : L 2 (R d ) ® J 7 -> L 2 (M d x Q, pgd^a; ® d/z) 

is unitary as well and satisfies U(ipo ® fi) = 1. A proof of the following 
lemma can be found in [8]. For completeness we sketch the argument. 

Lemma 2.4. T/ie operator Ue~ TLx U~ 1 is positivity preserving for T > 0. 

Proof. Define the cutoff function xn(u) '■= UX\y\<N- Abbreviate <\> '■= 4>{px) 
and set 4>n '■= 4>X\<p\<N- Since 4> is Hf bounded, one can show using the 
spectral theorem that Lq + Xc^n converges in the limit N — > oo in strong 
resolvent sense to L\. On the other hand Ue~ tLo U~ l is positivity preserving 
for t > (see XIII. 12 Example 3 in |32| or §1.4. in |34|). It now follows using 
the continuity of the functional calculus and the Trotter product formula, 
that 



Ue~ TLx U' 



lim 

N-too 



lim 

m— >oo 



jj e -TL /mjj-ljj e -TX<p N /mjj- 



Since Ue tx ^ N JJ 1 = e tx XN(£(px)) \ s a multiplication operator which is pos- 
itivity preserving, the lemma follows. □ 



8 



Theorem l2.3l holds as a simple consequence of Lemma l2.4l and Lemma [2,51 
below. 

Lemma 2.5. Let H be a self-adjoint operator in the Hilbert space L 2 (M, d/i), 
where fi is a probability measure. If e~ TH is positivity preserving, then 

lim --Log(l,e _rH l) = info-(iJ). 

T— >-oo 1 

Proof. The following proof is from [26]. Using that e~ TH is positivity pre- 
serving it follows that for / £ L 2 (M,dfi) with 

< ci < / < c 2 < oo, (2.3) 

the following inequalities hold 

cj(l,e- TH l)<(f,e- TH f)<cl(l,e- TH l). 

By Lemma 12.11 one concludes that E\ := inf supp [i\ = inf supp [if. Let X 
denote the set of functions of the form (I2.3P . The linear span of X is dense in 
L 2 (M, d/i). It follows that E\ = inf o~(H), since otherwise E\ > inf <r(-ff) and 
X(-oo,Ei)(H)L 2 (M,dfj,) would contain a nonzero vector which is orthogonal 
to X '. □ 



3 Path Integral Representation 

In this section we give path integral representations for 

Z T (P) :={n,e- TH ^Q), (3.1) 
Z T := (vo®tt,e- TLx ip ®n). (3.2) 

The path integral representation for (|3.1|) and (j3.2[) will be given in Theo- 
rem [XT] and Theorem 13. 2\ respectively. 

We use the following construction of measures, see [33] • We define the 
real Hilbert space N as the completion of 5 rea i(lR rf+1 ) with respect to the 
inner product 



{f )N = 2 r r f^,kmk ,k) ddkdko 

By the Minlos theorem, there exists a measure dP on 5^ cal (lR d+1 ) and a 
Gaussian random process £(/) on cS; cal (M d+1 ), which is indexed by / G N, 
with mean zero and covariance 

E rf p(C(/)e(5)) = (/,5)iV. 
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For / £ i?_i/ 2 and t £ R, we have f®6 t eN, and we set &(/) = £(/ ® <J t ). 
is a Gaussian random process, with mean zero and covariance 



VdpiUfMg)) = / f(k)g(k) 



2u)(k) 



e -u(k)\t-s\ d d k = (O,0(/) e -| s -*l^/ ( /,( ff )^). 



For f|3. If) . we introduce by (bt)t>o = {bi,t, bd,t)t>o Brownian motion 
starting at zero, i.e., a Gaussian process with 

^db{b a ,tbf3, s ) = S a) fiinhi(s,t). 

The following theorem is shown in [38] . For the convenience of the reader, 
we sketch a proof in the Appendix. 



(3.3) 



Theorem 3.1. For (|3.ip the following equation holds 
Zt{P) = ^db®dP 



exp(-A / i s (p hs )ds + iP ■ b T ) 
Jo 



For (|3.2p . we introduce qt = (qi,t, ■ ■, <ld,t) the oscillator process, i.e., a 
Gaussian process with 

1 



^dp(q a ,tq/3,. 



-5 a ^exp(-\t - s\). 



The following theorem can be shown in the same way as the Feynman-Kac 
formula was proven in [35], see also \26\ [9] and references therein. For the 
convenience of the reader we sketch a proof in the Appendix. 



Theorem 3.2. For (|3.2[) the following equation holds 

rT 



Zt — ~^dp®dP 



exp(-A I £t(p qt )dt) 
Jo 



(3.4) 



Remark 3.3. Using Kolmogorov's continuity theorem we can assume that 
t —7- £,t(Pbt) an d t Ct(Pqt) are a - e - continuous, and hence the integrals in 
(|3.4p and (|3.3[) exist a.e. as Riemann integrals. 



4 Power Series Expansion 

In this section we analyze the expansion of the r.h.s. of (|3.3p and (|3.4p . 
We use the BKAR decoupling formula to calculate the logarithm of these 
expansions. Our estimates on the resulting expression will allow us to take 
the limits (jl.4p and (|2.2p . In Subsection 14.11 we analyze the expansion for 
the harmonic oscillator. This subsection serves as a preparation for Subsec- 
tion H21 where the translationally invariant Nelson model will be studied. 
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4.1 The Expansion for the Harmonic Oscillator 

Integrating out the field in (|3.4p . we obtain 



Z T = E dp exp ( h& dP 



T i 2N 

A / &.(p qt )dt 



E 



dp 





T r T 



exp^A 2 ^ J W(q s -qt,s-t)dsdtj 



where we introduced the notation 

1 f\p(k)\ 2 



W{q,t) :=- 



\k\ 



e ^-i e -\M\t\ d d L 



(4.1) 



Expanding the exponential in a Taylor series, one obtains with g := A 2 /4, 



exp (A 2 / I W(q s - q t ,s - t)dsdt 
\ Jo Jo / 

\2n r rT rT rT rT 

1 + / / ••• / / W^-qt^si-h)... 

^ nl [Jo Jo Jo Jo 

x W(q Sn - q tn ,s n - t n )dsidti ■ ■ ■ ds n dt n 



Note that (14. lh is bounded, which allows to interchange integration and 
summation by dominated convergence. Inserting (|4.ip one finds, after inte- 
gration over the oscillator process 



oo n „ n n r n n / 

n^n^n^n (•■ 

n=0 J l-U> 1 \ j=l j=l JK j=l j=l v 

1 - 

— 'y ^ ki ■ kjAij(s,t) 



\kj\\sj tj | 



x exp 



(4.2) 



where 



and 



Aij(s, t) — C(si, Sj) — C(si, tj) — C(ti, Sj) + C(ti, tj) 



C(u,v)= l -e-\ u ~ v \ 

is the covariance of the one-dimensional oscillator or Ornstein-Uhlenbeck 
process. Clearly for any collections of times s and t, the matrix A(s,t) is 
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symmetric and positive-semidefinite. At this point one has achieved the 
mapping of the original model to a model of statistical mechanics. Indeed, 
Zt can be viewed as the grand canonical partition function of a gas of in- 
tervals on the real line (of course the bounds should be reversed if 
tj < Sj). These intervals are constrained by a two-body interaction corre- 
sponding to the ki ■ kjAij(s,t) terms in the exponential. Our first task is 
to write the quantity in (|4.2I) as an exponential. This is what cluster ex- 
pansions in statistical mechanics are typically used for. What follows is an 
adaptation of this general approach to the model at hand, using the BKAR 
formula. 

We will use the notation [n] for the set {1, 2, . . . , n}. For a finite nonempty 
set E, let us denote by E^ the set of unordered pairs {i,j}, where i and 



j are distinct elements of E. We will consider the space 



of multiplets 



u = [ui )i & e(2) indexed by I 6 E^ . A graph with vertex set E can be viewed 
as a subset of the complete graph E^ 2 ' . A tree on E is a graph without 
cycles which connects E. More generally, a graph without cycles is called 
a forest since it is a disjoint collection of trees, some possibly trivial, i.e., 
reduced to a single vertex. Let $ be a forest on E, and let h = (/ij)zg# be 
a vector of real parameters indexed by the edges I in the forest To such 
data one associates the element u($, h) = (u($, h)i)i €E (2) in E B(2) as follows. 
Let a and b be two distinct elements in E. If a and b belong to two distinct 
connected components of the forest 5, then u(J, M{a,6} = 0- Otherwise let, 
by definition, u($, h)s^ a ^ = mini hi where I belongs to the unique simple 
path in the forest J joining a to b. For E = [n] we define the function 



F(u) := exp 



^ n \ n 

-- ^2 kfAa(s, t) — - ^2 &i " kj Aij(s, t) x ttjjj} 



i=l 



»,J=1 



on the set M st2) . Observe that 



F(l) = exp 



1 n 

- ^2 ^» ' Aij{s, t) 



(4.3) 



where by 1 we denoted the point in 



with all components equal to 



1. Using the BKAR decoupling/interpolation formula [1 H [3] (see pQ for a 
pedagogical introduction) we have 

F(l) = [ dh (u($,h)) , (4.4) 

« f„™* ->M0,iP III 



5 forest ^I ' 1 ^ 
on E 



dui 
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where the sum is over all forests $ with vertex set E, the notation dh is for 
the Lebesgue measure on the set of parameters [0, 1]^ (here $ is identified 
with its set of edges), the partial derivatives of F are with respect to the 
entries indexed by the pairs belonging to 3, and the evaluation of these 
derivatives is at the h dependent interpolation points u(5, h). 

Note that each derivative 75-^ — , for {i, j} G produces a factor 

OU {i,j} 

ki ■ kj Aij(s,t). 
Besides, the content of the exponential now becomes 

■y n 1 n 

-^^2 k i A u(s,t) --^h-kj A i:j (s,t)u($,h){ id} . 

i=l i,j=l 

An essential feature of the BKAR decoupling formula is that it preserves 
positivity. This translates into the following lemma. 

Lemma 4.1. For all s,t G M. n , k G R dn , forest $ on [n], and h G [0, 1]% , we 

have 

-y n 1 n 

i=l *J = 1 

Proof. Since A a b(s,t) is positive-semidefinite, we have for any subset O of 

H 

h ■ kj Aij(s,t) > . 

i,j£<D 

For any given $ and h one can find (see PQ) nonnegative numbers Ai, X p 
satisfying Ylq=i A g = 1 as well as partitions n±, ...,7T P of [n], such that 

p 

q=l 

where for a partition ir of [n] we have defined (v^i := i{30 G 7r, Z C O}. 
Henceforth we will use !{•••} for the sharp characteristic function of the 
condition between braces. Now, we obtain 

Y n 1 n 
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1 n 

»=1 {mKN (2) 
^ n P 

9 fc ?^-ii( s ' *) + A( ? ( v ir q ){i,j} k i- kj Aij(s,t) 



i=l 



« =1 {i,j}£[n]W 



- ^ fc?A<(s, t) + Y, x iYl Yl • k 3 Avis* l ) 

q=l OGtT, {i,j}eO( 2 ) 



i=l 
P 



1 

oJ2 x « Yl S fci-M<j(M)>o 



□ 



We now insert (|4.4j) into 114. 2 j) and use Fubini's Theorem to pull the sum 
over the forest 5 and the integral over the parameters h out of the integrals 
over the s, t and k's. This is easily justified using Lemma 14. 11 

We therefore have 



n=l 



7t /I /i » ^ 

n n *i / n ^ 

5=1 i=i R " d J'=i 



i=i 



on [n] 
k 3 \\ Sj -t 3 \ \p( k j)\ 2 



\kj\ 



1 X II • % A/( s >*)) 



x exp 



j n j n 

o X] k i A ii( S i f )~ nJ2 ki ' k 3 Ai i( S > ^1 



i=l 



»J=1 



For given n > 1, we have a sum of the form which can be 

organized according to connected components: 

/(#) = — 1{ the J; form a partition of [n] } 



8 forest 
on In] 



k>l J u ...,J k C[n] 

E /(u* 

\{=1 

%l tree on J; 
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where the trees are summed over independently of each other. 

Inserting the above identity and noting the factorization which follows 
from the crucial property u($, h)^ j\ = if i,j belong to different compo- 
nents Ji, we have 



°° g n 1 
^ =1 + E^ElT E 



n\ *■ — ' k\ 

n=l k>X pi,...,pi->l 
Ei=iPJ=" 



1{VZ, I Jj| = pi} 11 { the ^ form a partition of [n]} 

Ji,...,JfcC[n] 
fe C 



«n e 

Z=l L X; tree on J; 



[0,l] s « 



1 n dt J 



{M. d )Pi 



x exp 



|fcj||Sj— tj\ 



lp(%)l 2 



(— fei - Ay Aij(s,t)) 



n An ( e 

1 1 n 



{»,3'}6Xl 



ij'GJi 



Note that we also used the fact u(^,h)^ jj = u(%,hi)^jy for a pair {i,j} 
in J/. Indeed, it follows from the definitions that the computation of this 
entry of the interpolation point is purely local to the component J/. 

It is easy to see that the contribution of a component «// only depends on 
its cardinality p\. Therefore, using the Multinomial Theorem and a trivial 
relabeling of the summation/integration variables, on can write 



00 1 1 
^ =1 + E^ElT E 



nl 



n\ ' k\ ' pi! ••■pi.! 

n=l k>l Pl,...,p fe >l 

T,LiPi= n 



n «" e 



X tree 
on \pi] 

Pi 



dh 



[o,i]i 



Pi Pz 
T ds 7 - TT dtj 



/ n a- n 

Pi 



g 1 1 1 s j % 



x exp 



Pi j Pi 

9 E^ 2j4 "( s '^ ~ 9 E fci ■ ^ ^j( s '*) u ( 3: '^){m} 



i=l 



#3 
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By cancelling the n\ and exchanging the order of the sums over n and k we 
immediately see that becomes the exponential of 



"**■==£*•£ 



n=l X tree 
on [nl 



dh 



n 

3=1 



x exp 



[o,ip 

|p(fc3)| 2 



[0,T] 2 



n n „ n 

n ^ n *i / n a- 



{M'}6X 



- n 1 n 



J} 



8=1 



i,3=l 



(4.5) 



Note that this last step relies on justifying the exchange of order of summa- 
tion which follows from the finiteness of 



71=1 

n 

n 

3=1 



x exp 



X tree 
on [n] 



dh 



[o,i]- 



[0,T] : 



n n „ n 

II- 1 ii * x_ II- - 

3=1 3=1 ^ 3=1 



- 1 kj 1 1 s j tj | 



x Yl dh - kj A i:j (s,t)\) 

{«}6I 



^ 71 1 n 

"n^^uCM) ~ ^z^h-kj Aij(s,t)u(%h){ id } 



i=l 



i,3=l 

«7^3 



The required bound IV < oo will be shown further below. 

In order to examine the limit limj^oo — lo s^ ^ j g b^ter t Q use [— ^ , 2] n 
instead of [0, T] 2n for the domain of integration over times. Now dividing 
by T and taking this limit essentially amounts to fixing a time say s\ to be 
the origin and replace the domain of integration over times by M. 2n . Indeed, 
the expression for the log in (]4.5p can be rewritten 



oo „ 

log Z T = y^l dsf- dt n K n (si, ...,t n ) 

, J\-T- T ]2n 
n=l J y 2 ' 2 J 

where the integrands are translation invariant in time (the sum over % is 
now included in the K n 's). So write 

K n (si, ...,t n ) = K n (0, s 2 - si, . . . ,t n - si) 
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and change variables to these differences 



log z T =^2 



OO „ T_ 

2 



ds 2 ■■■dt n K n (0, 



n=l 

AT T 

X 1 2~ Sl - S2 ~ 2~ Sl1 



Now change variables to u = ^ to find 



log Zt 



n=l 



ds 2 • • if n (0,S 2 , ■■■,t n ) 



xl{-T[- + u) <s 2 <T 



Note that the characteristic function goes to 1 pointwise almost everywhere. 
So, by the Lebesgue dominated convergence theorem, the key to taking the 
T — > oo limit is a bound on 



OO ,. 

V] / ds 2 ■ ■ ■ dt n \K n (0, s 2 , ■ ■ ■ ,t r 



or, which is stronger, a bound on 

oo „ 

r = / ds 2 ---dt n K n (0, s 2 ,..., t n ) 

where K n is defined in the same way as K n except one puts absolute values 
on the factors —ki ■ kj Aij(s,t) and g is replaced by its modulus. 

Note that the previous manipulations using translation invariance and 
the change of variable u = ^ easily show that TV < T x T. Therefore, 
modulo the key estimate T < oo which is proved below, we have shown that 



lim ■ 

T-j-oo 



log Zt 



n 

3=1 



^ git ^ 

n=l X tree 
on [n] 



dh 



[0,1] 



n n ,. n 

n da, n &i / n a 



\kj\\ Sj tj 



|P(%)| 2 



x (-ki ■ kj A i:j (s,t)) 



x exp 



- n 1 " 

yz2 k i A ii(.s,t) - - k i ■ k j A ij(s,t)u(1,h) {iJ} 



i=l 



i,3=l 



(4.6) 
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where, by definition, we set s\ := 0. 



4.1.1 Proof of the Key Estimate T < oo 

We first bound the big exponential by one, using Lemma 14.11 



n=l 



where 

r(")(£) 



n n ,. n 

n ds i n dt i / n A ; 



n 



7=2 7=1 7=1 
t ,| l^)| 2 



X fcj • fcj Ajj (s, t) I , 



with si = by definition. Note that the h parameters have been integrated 
out. Then one uses 

| hi ■ kj Aij(s,t)\ 



^ ~ k% ' \ kj\ ( 6 



+ e 



- | S j t j 



. +e -\u-s J \ +e -\u-t ] \j (47) 

Now one expands these factors so instead of X one now has a sum over 
refined trees (X, s) which is the data made of a tree X together with the 
knowledge, s, of which exponential decay one chooses for each edge: ss, st, 
ts or tt, that is s 6 {ss, st, ts, ft} 2 *. See the figure below for how this looks 
like. 
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The hard lines correspond to the edges of the tree X whereas the squiggly 
lines are 'internal' to its vertices. We can now, for fixed momenta k, perform 
the integrals over the times s and t. This integration is organized along 
the refined tree starting from the leaves and progressing towards the root 
corresponding to si which is not integrated over but set equal to 0. For 
the refined tree in the picture a possible order of integration is £4, S4, Sq, 
t®, tj, 87, ss, ts, tg, sq, ts, S3, t§, S5, S2, £2, and finally t\. The hard lines 
give factors of 1 by integrating the chosen exponential decay ie~l'"l in (|4.T[) . 
The squiggly lines produce j^y factors. Thus the refined tree (T, s) yields a 

contribution to r( n )(X) which is 

where di is the degree of vertex i in the tree T, as results from keeping track 
of the I ki I factors from (|4.7h . Then, summing over the choices of refinements 
s for a given tree %, it follows that 

r^(T) < 4 n " 1 JJ (^j 2\p(k)\ 2 \k\ d J- 2 d d k^j 

where we used that a tree with n vertices has n— 1 edges. Using the definition 
(jl.9p and for n > 2 we obtain 

r( n )(£) < 4™- 1 2™A 2n ~ 2 

where we used that d± + • • • + d n = 2n — 2. The bound in the n = 1 case 
only uses the hypothesis puj -1 £ L 2 . As a result, we find 

J^r (n) (T) < 4™-i2 n A n2 ~ 2 n n - 2 

by Cayley's Theorem for counting trees on n vertices. Then using n n ~ 2 < 
n\ e n ~ 2 , we finally have the desired estimate T < 00 provided |A| < Ao, 
where 

Ao = (2eA 2 )-5 . (4.8) 

4.2 Expansion for the Translation Invariant Nelson Model 

We integrate out the field in (|3.3p and we obtain 

Z T (P) = E db (exp(iPb T + A 2 jT W(6 S - 6 t , s - i))) • 
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We expand the second term in the exponential in a Taylor series and in- 
tegrate out the Brownian motion. Thereby we interchange integration and 
summation using dominated convergence and the boundedness of (|4.ip . As 
a result 



TP 

Z T (P) =e — — x 



£ 

n=0 



9_ 

111 



n n 



n ds i n *i / n a- 



3=1 



( e -\kj\\sj-ti\-{krP){sj-tj) 

3=1 



|P(%)| 2 



x exp 



1 n 

— ^ ^ ki ■ kj Aij(s,t) 



*,i=i 



(4.9) 



where g = ^, 

Aij(s, t) = C(si, Sj) — C(si, tj) — C(ti, Sj) + C(ti, tj] 

and 

C(n, t>) = min(u, v) 



(4.10) 



(4.11) 



is the covariance of the one-dimensional Brownian motion starting at the 
origin. Again, for any collections of times s and t, the matrix A(s, t) is 
symmetric and positive-semidefinite. 

We now follow the steps in §4.11 verbatim. The result can be summarized 
as follows. Using the same notations as in £ )4.1[ let 



< (\ki-kj Aij(s,t)\) 



n 

3=1 



e - 1 *i 1 1 »i -tj \-{k r P) ( Sj -tj ) \p( k j)\' 

\kj\ 



x exp 



^ n 1 n 



i=l 



*,3=1 
i+3 



where, by definition, we set s\ := 0. We now have the following result. 
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Proposition 4.2. Provided one can show the key estimate V < oo, we have 
that for any finite T, Zt{P) is the exponential of 



log Zt(-P) : 



n- 

3=1 



TP 2 ^ g» ^ /■ - 
2 ^ n! ^ A iii y 

n=l X tree" 7 ' ' 1 ! 
on [n] 



n n ^ / n ^ 

^ I j = \ j = l J ^ j = l 



\p(kj)\' 



x (-ki ■ kj Aij(s,t)) 



x exp 



-^^2 k i A u(s,t) - - y ki ■ kj A i:j (s,t)u(%h){ 



i=l 



»,i=i 

i¥=3 



Furthermore, under the same hypothesis, we have 
log Zt(P) P 2 A<? n 



lim 

r->-oo 



T 



on [n] 



72 1% n i v 

HdsjUdtj ]jd 



n 

i=i 



e _ | 1 1 Sj—tj \~(ky P)( Sj -tj ) \K k j)Y 



/ r„- „-i /— <t" 



x exp 



n 1 n 



(4.12) 



where, by definition, we set s± := 0. 0/ course, the A matrix now is given 
by Brownian motion instead of the oscillator process. 



4.2.1 Proof of the Key Estimate T < oo 



The beginning of the argument is the same as in §4.1.11 Indeed, Lemma [4.1 1 
holds for the new positive-semidefinite matrix A(s,t). Therefore, 



n=l V 7 X 
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where 



n lb lb n lb 

rW(x) : = / Y\ ds iY\ dt i T\d% 

x fl {e-toU'j-tihto-ni'i-tj) x H \-h ■ kj Ajis^l , 

with s± = and where the h parameters have been integrated out. 

Now is when we part ways with the argument of §4.1.11 We use, even in 
the case of a complex P, the crude bound 

e -{k r P){ Sj -tj) < e |fej||P||sj-^| (4-13) 

as well as the Cauchy-Schwarz inequality for the inner products of momenta 
in the factors attached to the edges of the tree T. This gives the inequality 



„ n n „ n 

r (n) (T)</ TldsjTldti Yld^j 



n 

x Y[ (e-NCH-PDK-M Ip^)! 2 !^-!^'" 1 ) x \Aij(s,t)\ (4.14) 

where, again, <ij is the degree of vertex i in the tree 1. However, we do not 
break |Ajj(s,t)| into four pieces as in f)4.T|) . Indeed the covariance C now 
no longer has exponential decay. The crux of our proof in the translation 
invariant case is that, nevertheless, the combination Aij(s,t) has a built- 
in decay which enforces the overlap of the time intervals for the indices i 
and j. This is a simple consequence of the independence of increments of 
standard Brownian motion. This makes the analysis similar to that for the 
spin-Boson model in [2]. For real numbers a, /3, let us use the notation 

I(a, (3) = [min(a, /?), max(a, /?)] (4-15) 

for the closed interval with endpoints a, f3, regardless of their relative order. 
We also use the notation |/| for the length of an interval /. Finally, let 
sgn(x) denote the sign of a real number x, or more precisely 

r 1 if x > , 
sgn(x) = < if x = , (4.16) 
( -1 if x < . 

We can now state the following crucial fact. 



22 



Lemma 4.3. For any s,t € W 1 , and i,j G [n], 

Aij(s,t) = sgn(si-tj) x sgn(sj -£,) x n I(sj,tj)| . 

Proof. One can check the formula using (|4.10p and (|4.1ip in all possible 
cases. Alternatively, one can notice that 

A ij (8,t)=E[(B a .-B ti ) (B Sj -B tj )} 

where B is the standard Brownian motion starting at the origin. Then one 
uses the independence of increments and the defining property 

E [(B u - B v f] =\u-v\. 

□ 

The main tool for the estimates is the following calculation. 

Lemma 4.4. For any real number fi > 0, integer p > and fixed real 
numbers a, (3, we have 

[ ds dt \s - t|P e -^-*l \I(a, (3) n I(s, t)\ = 2 (P + ~ "I . 

JUL* 

Proof. It is enough to treat the case where a = and (3 > which we now 
assume. Denote by I p the integral on the left-hand side. Let us first consider 
the case p = 0. By trivial symmetry 

J = 2 f ds dt e -^-*l |/(0, p) n J(s, t)\ l{s < t} . 

We now decompose the relevant integration domain given by the conditions 
s < t and 1(0, f3) n I(s, t) ^ into four pieces: 

I = 2(7 ,i + -fo.ll + Jo.lll + Jo.iv) 

where 



*>,I = / 

J 


ds dt e~^ 

R 2 


r-t| 


|J(0,/3) n J(s,t)| l{s < < f < /?} , 


Io,n= I 
J 


f ds dt e~^ 

R 2 


<-M 


\I(0,P)nI(s,t)\ l{0<s<t<(3}, 


Io,m = 

J 


f ds dt e - **' 1 

R 2 


r-t| 


|I(0,/3) n I(s,t)\ l{0<s<f3<t}, 


Io,iv = I 


f ds dt e - "' 1 


r-t| 


\l(0,(3)nl(s,t)\ l{s < and (3 < t] 



Jm. 2 
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Note that we ignored domain overlaps of Lebesgue measure 0. By 'time- 
reversal' symmetry we have Io,i = ^o,lll- Then by elementary calculus one 
obtains for these integrals the evaluations 

h,l = h,m = — 5-(e^-l-/x/9) , 
h,n = ^(2 + ^ + five"? - 2erf) , 



-'0,1V = 



Then taking the total we get 



M 2 



Finally, for p > 0, we use differentiation under the integral sign in order to 
show 



□ 



dfi ) Io ~{ dfi) ^~ vP+i 



We can now perform the integrations over times s and t, for fixed mo- 
menta k, in the right-hand side of (|4.14p . 

Lemma 4.5. We have, in terms of the vertex degrees d% in the tree X, 



/ 



n—1 



IP 'in ii' ' ' 1 x n \ A ^t)\ 

i=2 j=i j=i {ijyet 

3=1 J 

Proof. We use 1 6 [n] as the root of the tree X and orient the edges towards 
that root. The time variables are inductively integrated in pairs Sj, tj for 
each vertex j, starting from the leafs and progressing, following the edge 
orientations, towards the root. If j ^ 1 has vertex i as a parent, then the 
corresponding integral is 

f d Sj dtj e -IM(H-P|)h~t,l \ 8 . _ t^-i n I( Sj) tj)\ 
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= [\k,\(i - \p\)}^ 

by Lemma 14.41 The newly produced factor |sj — ti\ is dumped onto parent 
i. Since j has dj — 1 offsprings, the earlier integrations account for the 
\sj — featuring in the vertex j integral. The root requires a special 

treatment since s± = is not integrated over and the number of offsprings 
is d\ instead of d± — 1. In this case one has to compute 



[ dt x e-M(HP|)N \ tl 



<h 



2 di\ 



[|*i|(l-|P|)]*+i 

by a trivial Gamma Function evaluation. □ 

We now pick up the thread from Eq. (|4.14p and write, using the last 
lemma, 



r^(T) < 2 n (i - \p\y 3n + 2 x f\dj\ x ( ! 



d d k \mr 



By Cayley's Theorem which counts trees with fixed vertex degrees one has, 
for n > 2, 

n , s. n 

Enw- e (dl 4;.:. 

1 j=l d u ...,d n >l y ' y 1 j=l 

where we restored the T dependence in the notations. Therefore by the 
arithmetic versus geometric mean inequality 



d\ H Yd n 



X j=l d 1 ,...,d n >l 

£<l=2n-2 



<(n-2)!x ^ 

Ed;=2n 

<(n-2)!x2 n x ^ 1 

Ed 

< (n-2)! x 2 3n ~ 3 . 



di,...,d n >l 
Edi=2n-2 
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As a result, 

Therefore, the key estimate T < oo holds as soon as |A| < Ao(-P) with 



5 Convergent expansion for the effective mass 



We assume in this section that the cut-off function p(k) is rotationally in- 
variant, i.e., only depends on |fc|. From the formula ()4.12|) for the minimal 
energy E\(P) it is easy to obtain a double expansion in A and P. All one 
needs to do is expand all the factors e - ^'^^ - *^ as 



oo 1 



r,=0 r 



The previous estimates, and in particular the crude bound (I4.13p . show that 
this double expansion is convergent in the domain |A| < Xq(P). The degree 
zero term in P is of course E\(Q). It is also easy to see that the linear term 
vanishes. Indeed, such a term contains integrals of the form 

„ n 

/ \[d% [-{k h -P){s h -t h )]xW{k) 
Jm. nd - =1 



where 



w{k) = [] ( e -i*-M 



\p(kj)\ 2 



x (-ki ■ kj Aij(s,t)) 



x exp 



-y n 1 Tl 

■^^2 k i A u(s,t) - -^ki-kj A ij (s,t)u(%h)^ id} 



i=l 



»,J'=1 
i+3 
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and where we suppressed the dependence on n,T, s,t for lighter notation. 
By rotation invariance of p, one has that W(k) is invariant by simultane- 
ous rotation of all momenta kj. This implies the vanishing of the previous 
integral. As for the quadratic term in P, it is clearly equal to 



y y y 

n=1 jr'l,jr'2 = l ^l,l / 2 = l 1 tree 1,1 

on [n] 



dh 



i n ds J n dt i 

J=2 j=l 



« n 

/ T kj ltVl kj 2 ,v 2 Py 1 Pv 2 ( s ji ~ tjl)( s j2 ~ ^32) X W(k) . 

Using the invariance property, and a change of momentum variables which 
is a reflection in the z^i-th coordinate, it follows that only terms with v 2 = 
v\ = v survive. Besides, by rotation invariance, these are the same for all 
v. Therefore, the quadratic part is of the form 2 m s (A) w ith 



1 OO \ 2/1 ^ f n lb II: 

-4 = i- Fr^i T T dh ir,/s ir,// 

m e g(X) d^^nl ^ ^ L ir i Lan-i 11 J 11 ' 

cttv y n=l jij'2=lX tree" 7 ' ' 1 ! JK j=2 j=l 

on [n] 

„ n 

/ I ^ (%l " ^'2) — tjl)\ s jl ~ ^2) 

iR" d 7=1 



i=i 



g 1 1 s j *J° I 



\Kkj)\ 2 



x (-fci • fej Aij(s,t)) 

{i,j}ei 



x exp 



^y.kfAuis^) - ~y^h - kj A i:j (s,t)u(%h) { 



i=l 



*J=1 



One can then trivially substitute into the geometric series = 1 + x + 
x 1 + • • • in order to obtain a convergent expansion for the effective mass 
m e fr(A) itself. 

It is important to remark that the first nontrivial term with n = 1 is 
nonvoid, since there exists a tree T on the singleton {1}: the empty tree! 
Indeed, if one writes the expansion for the inverse mass as 



1 

m e ff(A) 



1 + c 2 X 2 + c 4 A 4 + • • • 
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then one can easily extract the coefficient C2 as follows. In this case one has 
n = 1, ji = j2 = 1, and T = 0. There is no integral over h parameters and 
our formula reduces to 



f-2 



1 

u 



dU 



d d h \h\ 2 tje 



|*i||ii| |p( fc i)l ; 



exp 



1 



which can readily be simplified to 



C2 



d d k 



\p(k)f 



1 + 1M 

1 ~ 2 



Therefore the expansion for the mass is m c g- (A) = 1 — C2A 2 + 0(A 4 ) with 
positive A 2 coefficient. Indeed, the physical picture is that the particle gets 
heavier when it is dressed by the Boson cloud. 
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Appendix 

In this appendix we sketch proofs of Theorems 13.21 and 13.11 Without loss of 
generality we can assume that A > 0. First we address Theorem 13.21 Define 
the cutoff function xn{u) = max(y, —N). Set 

4>x,N = Xn(4>(Px))- 

Since 4>(p x ) is Hf bounded it follows by the spectral theorem that := 
H osc + Hf + X(px,N converges in strong resolvent sense to L\. Hence 

(ipo <g> n, e~ TLx ip <g> tt) = lim (ip <g> tt, e~ TLx - N ip ® fi). 

7V->oo 
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Using the Trotter product formula, we obtain 

(ip ® n, e- TL ^ N ip ® fi) = lim (c^ ® O, ( e - ri o/m e -TA^^/ m | m ^ 3 ^ 



lim 

m— >oo 



exp(-^A X 7v(^(p< ?tl ))Ai) 



i=l 



E, 



dp®dP 



T 



exp(-A / XN{it{p qt ))dt) 
o 



where in the second equality we used the FKN formula, see |34| . with the 
notation At = T/m and tj = iAt. In the third equality we used dominated 
convergence and that the sum converges to the integral a.e., since t h-> £t(Pq t ) 
is a.e. continuous. Now we use monotone convergence, to conclude that the 
right hand side converges to 



Erf p ®rfp(exp(-A / £t(pq t )dt)), 
Jo 



as N — > oo. This shows Theorem! 

Now we sketch a proof of Theorem 13.11 In contrast to [38] , the proof 
below uses a Gaussian measure to linearize a square and also the so-called 
pull-through formula. Using an analogous argument as above we have 

(n,e- TH ^n)= lim (ft, e -™( p )+ A <^)ft). 

N— s-oo 

Using the Trotter product formula, we obtain 

(fi,e- T( " o(P)+M) '^) 



lim (n, \ e -TH {P)/m e -TX4> 0:N /m\ m Q) 
m— >oo I J 



E 



lim f dp m (a)(n,f\ \ e ^AP-Pf) e - TH f/ m e- TX ^/ m \n) 

J 3=1 

/m 
dp m (a)(U, e^ p TT \e- TH fl m e- TX4, ^ N/m X fi) 
- LJ - I J 

j'=i 

m 

e i6 - p exp(-A^Xiv(^(p 6% ))Ai) 

e i6r F exp(-A [ T XN(UPb t ))dt) 
Jo 



lim E, 



db®dP 



db®dP 
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where in the second equality we introduced the Gaussian measure dfi m on 
]R 3m B a = («!, ...,a m ) with mean zero and covariance T/m. In the third 
equality we used the pull-through formula and the notation f3 n = Y^j=i a j- 
In the fourth equality we used again the FKN formula. The last equality 
follows from the a.e. continuity of paths and dominated convergence. Now 
the right hand side converges to 

E d 6®dp(e* r-P exp(-A f Zt(p bt )dt)) 
Jo 

as N — > oo. For P = this follows from monotone convergence, then it 
follows for P 7^ by dominated convergence. This shows Theorem 13.11 
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